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ABSTRACT. Due to the fact that ZIKA virus circulate in both primates and human beings 
populations, the dynamics of the virus can be described better by a mathematical model 
which include both populations. This work presents a compartmental mathematical model, 
SEIR-SEI-SEIR, for human beings, mosquitoes, and primates respectively. ‘The model con- 
siders sexual transmission for both human beings and primates populations, and vertical 
transmission for all the three populations. The work examines effects of some model pa- 


rameters on basic reproductive number, on human beings disease incidence, and on human 





beings disease prevalence. Our study finally concludes that neither vaccination of human 
beings alone nor the application of mosquito nets, condoms, and abstinence from sexual 
activities alone can eradicate the disease. Instead, one can eradicate the disease by apply- 
ing one, or a combination, of the following possible, but not limited to, sets of strategies: 
vaccine both primates and human beings; heavily use of mosquito pesticides that can let 
mosquito’s death rate be at least 15 times more than its natural death rate; a combination of 
the use of mosquito pesticides, and the use of mosquito nets, condoms, or abstinence from 
sexual activities; a combination of the use of mosquito pesticide and vaccinating human 


beings /primates. 


1. INTRODUCTION 


ZIKA virus is a mosquito - transmitted flavivirus of the family fraviviridae [13]. It causes 
fever, Guillain-Barre syndrome (GBS), and Microcephaly - brain defects to a fetus of an 
infected mother [12]. The congenital neurological disorder due to ZIKA disease poses a 
significant concern for the nearly 15 million children born in the Americas each year|17]. 
Moreover, WHO estimated that about four million Americans could be infected by ZIKA in 
2016 [10]. 

ZIKA virus is primarily transmitted through mosquito bites of an Aedes genus |11]. It can 


also be transmitted through sexual intercourse |9, 15]. Vertical transmission from an infected 





mother to her fetus during and before delivery is possible too [1]. Infected mosquitoes can 
also lay eggs which are infected by ZIKA virus [18]. 
ZIKA virus was first isolated from sentinel rhesus in 1947 in Uganda at Zika forest [14]. In 


Nigeria, the first human infection was reported in 1954; thereafter, ZIKA outbreak occurred 
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in Yap Island in the Federated States of Micronesia in 2007 [14]. This was followed by a 
large outbreak in French Polynesia in 2013-2014, then it spread to New Caledonia, the Cook 
Islands and Eastern Islands and then to Brazil in 2015 [14]. 

Many studies on ZIKA dynamics have been conducted since the 2015 ZIKA outbreak in 
Brazil. However, we lack a comprehensive mathematical model which consider dynamics 
of ZIKA virus among and between human beings, mosquitoes, and primates populations. 
This work presents a compartmental mathematical model with human beings and primates 
as the virus hosts, and mosquitoes as vectors. The model includes sexual transmission for 
both human beings and primates populations, and vertical transmission for all the three 


populations. 


2. MATERIALS AND METHODS 


2.1. Model Development. We regard human beings and primates as the current known 


hosts of the virus, while mosquitoes (Aedes species) act as transmission vectors. We use the 





SEIR, SEI, and SEIR compartments for human beings, mosquitoes and primates respectively 
to model the dynamics of ZIKA virus disease. We neglect the human death due to ZIKA 
since very few cases have been reported and not yet confirmed by serious examination. We 
assume that human beings and primates are the only hosts of ZIKA virus. We also assume 
that all Aedes mosquitoes responsible for ZIKA transmission have similar transmission prop- 
erties, and therefore, we categorize them in one population. We assume that the epidemic 
process is deterministic. We ignore migrations in all the three populations by assuming that 
emigrations counter balance immigrations. 

In our model, a susceptible human (Sn) gets the infection in two ways: first, through 
having sex with an infected human(J;,) and second, through a bite of an infected mosquito 
(Im). He/she then progresses to an exposed class at a rate which is proportional to contact 
rate and level of infections. However, not all contact can lead to transmission of the infec- 
tion. Transmission of infection depends on susceptibility of the susceptible, infectivity of the 
infected member, and probability that a contact will transmit infection [5]. 

In this model, susceptible members leave the susceptible class to exposed class at a rate 
By (Anntn/Nit+tAmntm/ Ne), where, 3; is human being susceptibility, Ap, is an adequate human 
to human contact rate, and Amp is an adequate contact rate from mosquito to human. 
The fractions [j,/N,, or Im/N1ı, represents a possibility that a person or a mosquito that is 
contacted with a susceptible is infectious [5]. 

Since ZIKA virus can be transmitted vertically, there is a proportional pı that a new born 
in human population N; is born susceptible at a rate y,; and a proportional 1 — pı that the 
new born is infected. A susceptible human beings remains in an exposed class for some days 
after he/she gets the virus. An exposed human (£;,) progresses to the infectious class at a 


rate a; and then recovers naturally at a rate ag and attains permanent immunity against 
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reinfections. Since the death due to ZIKA disease is very rarely, almost negligible, the model 
considers only natural death at a rate uı for human. 

A susceptible primate (S,) gets infection in two ways: first, through having sex with an 
infected primate (Jp), and second, through a bite of an infected mosquito (Im). It then pro- 
gresses to an exposed class at the rate $3(Amplm/No2+Applp/N3). Due to vertical transmission, 
there is a proportional p that a new born in primate population N3 is born susceptible at a 
rate y3, and a proportional 1 — p that the new born is infected. An exposed primate (Ep) 
progresses from an exposed class to infectious class at a rate a3. It then recovers naturally at 
arate as; and attains permanent immunity against reinfection. A primate dies naturally at 


arate u3. No death of a primate due to ZIKA virus has been reported. A susceptible Aedes 
TABLE 1. Model Variables 


Susceptible human beings 
Exposed Human beings 
Infected Human beings 
Recovered Human beings 
Susceptible Mosquitoes 
Exposed Mosquitoes 


Infected Mosquitoes 


Susceptible primates 
Exposed primates 
Infected primates 


Recovered primates 





~ ~ ~ 
pra pe harm & 





mosquito (Sm) gets the infection when it bites either an infected human (/;,) or an infected 
primate (/,). It then progresses to an exposed class at the rate B2(Apmln/Ni + Apmlp/N3).- 
A proportion, pə of eggs laid by an infected mosquito are not infected. Meanwhile, the re- 
maining proportion 1 — p> of eggs laid by an infected mosquito are infected. Approximately 
1 egg out of 290 eggs of infected mosquito are found to be infected with ZIKA virus [18]. 
Mosquitoes are recruited through birth at a rate yo. ‘The susceptible mosquito remains in 
an exposed class for some days and then progress to infectious class at a rate ag. It then 
dies naturally at a rate u3 after about 1 month from the day it hatched. Table 1 and Table 


2 show model’s variables and parameters respectively. 
Table 2: Model Parameters 
Parameter Description Units 


Ua Natural death rate of a human being death day™t 


Continued on next page 
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Table 2 — Continued from previous page 


Parameter Description Units 
Uo Natural death rate of a mosquito death day~! 
U3 Natural death rate of a primate death day~' 
v4 Birth rate of a human being birth day~* 
Vo Hatching rate of a mosquito hatches day~' 
"V3 Birth rate of a human primate birth day~* 
pı Probability that a human being is born infected unit less 
D2 Probability that a mosquito is hatched out infected unit less 
D3 Probability that a primate is born infected unit less 
By Human being’s susceptibility unit less 
Bo Mosquito’s susceptibility unit less 
Bs Primate’s susceptibility unit less 
Oy Human being’s progressive rate from exposed class to 

infectious class humans day~! 
Q92 Human being’s recover rate humans day~* 
Q3 Mosquito’s progressive rate from exposed class to 

infectious class mosquitoes day! 
a4 Primate’s progressive rate from exposed class to 

infectious class primates day™! 
Os, Primate’s recover rate primates day! 
Ahh Adequate contact rate that lead to infection from 

human being to human being through sex contacts day~* 
App Adequate contact rate that lead to infection from 

primate to primate through sex contacts day~! 
Ahm Adequate contact rate that lead to infection from 

human to mosquito contacts day~' 
Arah Adequate contact rate that lead to infection from 

mosquito to human contacts day~* 
Apm Adequate contact rate that lead to infection 

from primate to mosquito contacts day~' 
Arip Adequate contact rate that lead to infection from 

mosquito to primate contacts day~* 





The dynamics of ZIKA described in Section 2.1 can be represented by the compartmental 
diagram shown in Figure 1. The doted line shows interactions that lead into transmission 


of the infection, and the arrows shows movements of individuals, mosquitoes, and primates 
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from one compartment to another. 








HiSh, 


A In 
By (Ann N. + Amn 


P1ViN, 


(1 — p2)y2N2 





(1 — p3)y3N3 


FIGURE 1. Flow diagram of ZIKA infections between human beings, 


mosquitoes, and primates. 


2.1.1. Model Equations. The model equations as deduced from the flow diagram in Figure 
1. The equations are categorized into three groups: human beings, mosquitoes, and primates 


populations. 
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Human beings 


dS 
(1) aE = pi — Se — BiAnntn/N1 + Amida Na) Oh 
dE}, 
(2) Pas Pal Arnd Nit Amntm/No) Sh — HEr = 1 En, 
dl, 
(3) TE = aL, + (1 E Pym — HWıln — Q&2lh, 
dR 
(4) ara = QoL), = ui Rp. 
Mosquitoes 
Gielen 
(5) Po = pyy2N2 — beoSm — B2(Anmta/M1 ac Anmol No lOni 
dEm 
(6) dt — Bo(Anmtn/™1 + Ned N3) Om E Ho Em z a3 Em, 
dlm 
(7) P az Em + (1 — p2) Y2 N2 — Holm. 
Primates 
AS 
(8) dt = P33 N3 _ [39 _ B3(Amplm/N2 a Applp/ N3) Sp, 
dE 
(9) a5 = B3(Amplim/Na + Applp/N3)Sp — U3 Ep — a4 Ep, 
dI 
(10) A = Q4Ep + (1 — p3) Ys Ns — Hlp — @sIp, 
dR 
(11) a = asl, = l3 Pp, 


wheres Ny =r t- Det Iet Rig N = Sayt Bala N SS Sp A Rad 
CRE EE ey ae E a) 


3. MODEL ANALYSIS 


3.1. Invariant Region. Our system of differential equations for human being, mosquito 
and primate, all together can be written in the form dX/dt = Aq) X + F. 

Ko (Sis ig oma El es Ops peli) 3 

F = (p11, 0, (1 — pr) 71N1, 0, poy2No, 0, (1 — pe) y2No, pays N3, 0, (1 — p3)ys N3, 0)", and 
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—a 0 0 0 © 0 0 0 0 0 0 
iG b 0 0 0 0 0 0 0 0 0 
0 a —c 0 0 0 0 0 0 0 0 
0 0 œa -m 0 0 0 0 0 O0 0 
0 0 0 0 -d 0 0 0 0 0 0 

(12) Awy=} 0 0 0 0 kb -e 0 0 0 0 0 |, 
0 0 0 0 0 œ -m 0 0 O0 0 
0 0 0 0 0 0 0 -f 0 0 0 
0o 0 0 0 0 0 0 Bago 0 
0o 0 0 0 0 0 0 0 œ -h 0 
co 0 © 0 0 0 0 0 a =p 


where: a = fy + kıĝı,b = fy + @1,C = pa + O2,d = u2 + koßb2,€ = po + Q&3, f = u3 + 
k383, g = u3 + a4, h = u3 + as, kı = Anntn/ Ni + Amntm/ Ne, k2 = Anmtn/N1 + Apmtp/N3, and 
k3 = Amplm | N2 + Applp/ N3. 

Since A(x) is a Metzler matrix, i.e. a matrix such that off diagonal terms are nonnegative for 
all X, and F > 0, we conclude that the system of differential equations for human beings, 


mosquitoes, and primates is positive invariant in Ry [8]. 


3.2. Positivity of Solution. So long as we are dealing with populations, all variables 





(solutions) in the system of the differential equations must be greater or equal to zero (non- 
negative) for t > 0. We first let the initial values (i.e. at t = 0) of the variables as follows: 
(En(0), Em(0); Ep(0), Ln(0),Lm(0), Lp(0), Rn (0), Rp(0)) > 0, (Si(0), Sm (0), $p(0)) > 0. Con- 
sider equation (1): ar = py yi Ny — HSn — Bi(Anntn/ Ni + Amh Im/N2) Sh. Then, 

= > =M Sn — B1(Analn/Ni + Amntm/No)Sy. Solving the differential inequality we get a 
solution, S;,(t) > Ce Wit 81Annta/Nit+BirAmntm/N2a)t Put initial conditions: t = 0 and S(t) = 
S,(0), we get C = S,(0). Thus, we get S(t) > Sp_(O)eT Mt tAnnta/Ni+BirAmntm/Na)t Since 
ea +Birnnth/Ni+FirAmntm/Na}t > Q, then, S,(t) > 0, Vt > 0. Using the same approach, we have 
Se) 0S) 0.5.) 0B) 0b 0. BG) 2 OL) SO.) S01.) 
0, R(t) > 0, and R,(t) > 0 for Vt > 0. 


3.3. Disease free equilibrium point. We find the disease free equilibrium points for 





three populations, human beings, mosquitoes, and primates, separately before combining 
the points to get the equilibrium point for the whole system. 


Consider system of differential equations (1) - (4) that represents the dynamics in human 
a5, = dEn — dy dRe 


beings population. At the equilibrium point = = = = Ue For disease free 


situation Ep = Ip = Rp = 0, Sh = Ni, p = 1. Since Se = 0, then an = 0, thus N; = Ny = 
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constant. From equation (1), we have: 


(13) i pıyı NT / ea 


Substituting pı = 1 into equation (13), we get Sp = vii, Following the same procedures 














for mosquitoes and primates populations, we obtain Sm = zi and Sp = pi, Thus, the 
disease free equilibrium point for the whole system is (nS „0,0,0, a , 0, 0, rå; ,0, 0,0). 


3.4. Basic reproductive number. Basic reproductive number is the number of sec- 
ondary infections produced by an infectious, during his/her entire infectious period, when 


exposed to a purely susceptible population [16]. When the basic reproduction number is 





greater than one, the disease-free steady state of a dynamics system is unstable, and when 
it is less than one the infection-free steady state is stable |6]. 

We use the next generation method to find the basic reproduction number of ZIKA dis- 
ease. Consider system of differential equations (1) - (11). Let f; be a vector that consists of 


new infectious terms of transmitting components and v; be a vector of the remaining transfer 


terms. 
Bi (Anntn/N1 + Amntm/Ne)Sh py Ep + iEn 
(1 — pi) N — 0 En + path + Aan 
0 —agl, + pp, 
Fis Bo An lng Nt F Apmp Na) Sm mE H2 Em + A3 Em 
(1 — pa) y2 N2 1 —03 Em + Halm 
B3(Amplim/N2 + Applp/ N3) Sp u3 Ep + a4 Ep 
(1 — p3)y%3 N3 —&4 Ep + Ugly + As5ty 
0 —Qs lp + u3 Pp 





Let F and V be the Jacobean of F; and V; respectively, evaluated at the disease free equi- 


librium point. 


0 Bàn 0 0 Ym 9 0 0 
0 0 0 0 0 0 0 0 
0 0 0 0 0 0 0 0 
Bor m N2 B2Apm N2 
aap? e C 0 0 apn’ 0 
0 0 0 0 0 0 0 0 |’ 
0 0 0 0 Bae 0 B5rpp 0 
0 0 0 0 0 0 0 0 
0 0 0 0 0 0 0 0 
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mitar © 0 0 0 0 0 0 

0; pitas 0 0 0 0 0 0 

0 zsa Mo O 0 0 0 0 

7 0 0 0 Aas Ü 0 0 0 
0 0 0 =i m 0 0 0 

0 0 0 0 0O u3+as 0 0 

0 0 0 0 0 —a, uz+a; 0 

0 0 0 0 0 0 =e ae 


The basic reproductive number is the largest Eigenvalue of a matrix FV~! and is given by: 


3C — B? 
(14) Ro = 1/6P — 2/3(——5—) - 1/3B. 
P = (—108D + 36C B + 12v12B3D — 3B?C? — 54BCD + 12C3 + 81D? — 8B3)1/3, 
= = = = Girne a= B3Appo 
Boa Ht pC = Op ee SM Ok aCe), 0 = ee ee P Gee ee 


BiAmnNia3 p= B2AnmN2a01 kn B2 Apm N204 m = 3 Amp N3Q3 
Nə(u2+03)(u2)’ Nı(uı+a1)(mı+a2)’ N3 (u3 +a4)(u3+as5)’ N2(u21+a3) ("2 ` 


C= 
3.0. Parameter estimation and simulations. Most of parameter values used in this 
study were obtained from different studies on ZIK A disease. Some parameters were estimated 
and some were assumed. The reader should also note that all parameters used in this study 
are subject to change. Mosquito biting rate for example, may vary according to the weather. 
Human-human sexual contact rate may also vary from one group of people to another. 

We consider the adequate contact rate as the product of a contact (or a biting) rate and 
a probability of an infected one to transmit infection. We estimate that a human being 
can make an average of 2 sexual contact per week (i.e. 2/7 per day) and primate does 
make an average of 1 sexual contact per week (i.e. 1/7 per day). We estimate a range, 
0.1 — 0.75 as a possibility of infected person to transmit infections, and we assume the 
same for infected primate. ‘Thus, we calculate adequate contact rate from human to human, 
Ann = 0.6(2/7) = 0.17 per day, and the adequate contact rate from primate to primate, 
App = 0.6(1/7) = 0.086 per day. 


Table 3: Model Parameters 


Parameter Range Value Sources 

[Ly 1/21900 [4] 

H yu A 

H3 1/5475 [2] 

yı 1/21900 Assumed = pı 
Y2 1/14 Assumed = [lo 
13 1/5475 Assumed = u3 


Continued on next page 
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Table 3 — Continued from previous page 


Parameter Range Value Sources 
pı 0.8-1 0.95 Estimated 
D2 0.997 [18] 

P3 0.8-1 0.97 Estimated 
By 0-1 0.85 Estimated 
Bo 0-1 0.96 Estimated 
D3 0-1 0.85 Estimated 
a1 1/7 - 1/2 1/4 [10] 

Q92 1/14 - 1/3 Le [3] 

Q3 1/12 - 1/8 1/8 [10] 

4 1 /6-1/2 1/4 [7] 

Os at most 1/9 1/12 [7] 

Anh (0.1 - 0.75) (2/7) 0.17 Estimated 
Ap (0.1 - 0.75)(1/7) 0.086 Estimated 
Ahm (0.1-0.75) (0.3-1) 0.18 [3] 

Arah (0.1-0.75) (0.3-1) 0.21 [3] 

Apm 0.30 [2] 

Ann 0.35 [2] 


The susceptibilities, 61, G2, and 83 (possibilities of getting infection after being exposed to 
the adequate contact rate) for human being, mosquito and primate are assumed to be 0.85, 
0.96, and 0.85 respectively. The assumptions base on the fact that susceptibility, herein, 
depends on natural immunity against ZIKA infections, and we roughly estimate that very 
few (15% of human beings population, 4% of mosquitoes population, 15% percent of primates 
population) can resist the infection naturally. Thus, 85%, 96%, and 85% are the respective 
possibilities for human being, mosquito, and primate to catch the infection, when exposed 
to the adequate contact rate. For the purpose of simulations, we assume that we deal with 
stable populations that birth rate is approximately equal to death rate as shown in ‘Table 3. 

Recovery rates, for both primates and human beings populations, are approximated to 
be equal to the reciprocal of their respective infectious periods, as it is proved to be an 
appropriate mathematical assumption in [5]. We hold the same assumption and let the 
progressive rates from the exposed classes to infectious classes be equal to the reciprocals of 
the latent periods. We roughly estimate the possibility of a person, and that of a primate, 
to be born with infection as 0.05 and 0.03 respectively. ‘This is to say, out of 100 live birth, 
5 babies are born with ZIKA infection, and 3 primates out of 100 live birth born with the 


infection. ‘Table 3 shows parameter’ ranges, values, and their origins. 
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Variation of Human beings population with Time Variation of Primates population with Time 
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600 Variation of Mosquitoes population with Time 
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FIGURE 2. Variation of populations with time. 


In order to understand the behaviors of variables of the model, we plot the graphs of 
the variables against time. Figure 2(a) shows the variation of human beings population (in 
four classes: susceptible, exposed, infectious, and recovery) with time. while the suscepti- 
ble human beings always decreases from initial value and approach zero asymptotically, the 
recovered human beings always increase and seems to approach a constant value asymp- 
totically. On the other hand, exposed human beings and infectious human beings initially 
increase from their starting values, attain their maximum values, and then gradually de- 
crease and approach zero asymptotically. This is so because we set the birth rate equal to 
death rate, and recovered individuals are not going back to susceptible class since they attain 


permanent immunity against new infections. Therefore, as time gets larger and larger the 





susceptibles decreases since new birth counter balance with death. The recovered individuals 
increase as susceptibles decreases because we regard ZIKA as not a fatal disease, and hence 


those who leave the susceptible class ultimately end to the recovered class. 
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Disease Incidence 
Disease Prevalence 





0 20 40 60 80 100 120 
Time[days] Time[days] 


(A) (B) 


FIGURE 3. ‘The plots of disease incidence and disease prevalence against time 


for the different values of susceptibility of human beings. 
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FIGURE 4. ‘The plots of disease incidence and disease prevalence against time 
for different values of human to human adequate contact rate through sexual 


activities. 


Figure 2(b) shows the variation of primates population (in four classes: susceptible, ex- 
posed, infectious, and recovered) with time. The variation of primates with time exhibits 
similar behavior as observed in human being population. 

Figure 2(c) shows the variation of mosquitoes (in three classes: susceptible, exposed, 
and infectious) with time. The infectious mosquitoes, and exposed mosquitoes, initially 
increase to its maximum value and gradually decrease and approach zero asymptotically. 
On the other hand, the susceptible mosquitoes decreases as the infectious mosquitoes in- 


creases, and increases as the number of infectious mosquitoes decreases. This is because, 
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FIGURE 5. Variation of Rọ with adequate contact rates. 


in mosquitoes populations there is no recovery class. Infectious mosquitoes will remain in- 
fectious throughout their entire life. Hence, all mosquitoes who are leaving the susceptible 
class will ultimately end to the infectious class. 
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The Graph of Basic reproductive number Vs Human Susceptability 
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FIGURE 6. Variation of R, 
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The Graph of Basic reproductive number Vs Primates recovery rate 
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FIGURE 7. Variation of Ro with recovery rates and mosquito birth rate. 


ZIKA VIRUS AMONG HUMAN BEINGS, PRIMATES, AND MOSQUITOES POPULATIONS 15 


It is also very important to understand the effects of varying some parameters on disease 
incidence, and on disease prevalence. Figure 3(a) shows disease incidence, in human beings 
population, decreases with a decrease of human susceptibility. Figure 3(b) shows disease 
prevalence, in human beings population, also decreases with a decrease of human suscepti- 
bility. This imply that number new infection, and proportion of people who are infected can 
be reduced by reducing human susceptibility. Remember, susceptibility can be reduced by 
vaccinating some proportion of members of a population as described earlier. For example, 
Figure 3(a) shows that when susceptibility, 8; = 0.85, the peak value of new infection rate is 
0.025 incidences per day, and when susceptibility, 6, = 0.25, the peak value of new infection 
rate is about 0.006 incidences per day. 


Figure 4(a) shows disease incidence, in human beings population, decreases with a decrease 





of human-human adequate contact rate through sex. Figure 4(b) shows disease prevalence, in 
human beings population, also decreases with a decrease of human-human adequate contact 
rate through sex. This suggests that the use of condoms or abstinence from sexual activities 
can help to reduce both new infection and proportion of people who are infected by ZIKA 
virus. 

To understand more about the impact of varying the parameters on existence of the 
disease, we plot graphs of basic reproductive number against parameters. Figures 5(a) 
- 5(f) show that the basic reproductive number increases as the adequate contact rates 
increase. For example, Figure 5(a) shows that reducing human adequate sexual contact rate 
(either by using condoms or abstinence from sexual activities) lead to a decrease of the basic 
reproduction number. However, the plot shows that when the adequate sexual contact rate 
is zero, the basic reproductive number is still greater than one. This is to say, reducing 
sexual contacts even up to zero does not result in eradication of the disease. 

Moreover, when we tried to vary parameters, by keeping all human adequate contact rates 
(human to human, human to mosquito, and mosquito to human) zero, and left non-human 
adequate contact rates and other parameters with their values shown in Table 3, the basic 
reproductive number was still greater than one. The observation suggests that both human 
beings and primates populations should be considered if we want to eradicate the disease 
by reducing adequate contact rate alone - through application of mosquito nets, condoms, 
and abstinence from sexual activities. But in practice, it is impossible to let primates use 
mosquito nets and condoms or abstained from sexual activities. Thus, we need to incorporate 
other control strategies such as vaccination and the use of pesticide that kills mosquitoes 
rather than relying on mosquito nets, condoms, and abstinence from sexual activities alone. 

It should be noted that the above conclusion is true only if the values App, Amp, Apm, 03, 
and a, are the same or greater than as they are estimated and appeared in Table 3. If 
we deal with a primate population in which these parameters are much less as appeared in 
Table 3, for example, we can possibly eradicate ZIKA disease through the application of 


mosquito nets, condoms, and abstinence from sexual activities alone. For example, if we 
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set App = 0.001, Amp = 0.001, Apm = 0.001, 63 = 0.71, and a4 = 0.001, and reduce App, Ama, 
and Aam to 0.01, 0.1, and 0.1 respectively, the basic reproductive number reduced to 0.7434 
which implies the disease disappearance. 

Figure 6 shows the variation of the basic reproductive number with susceptibilities. ‘The 
basic reproductive number increases as a susceptibility of human or primate increases. As 
we described earlier that susceptibility is a parameter that reflects the natural immunity 
against ZIKA infections one can possess i.e. a population with no individual who possesses 
natural immunity against the disease has 1.0 susceptibility, and a population in which 20% 
of its members possess the natural immunity has 0.8 susceptibility. 

Thus, susceptibility can be reduced by increasing natural immunity against ZIKA infection 
by possibly, vaccinating some proportion of either human beings population or primates 
population or both. By doing so the basic reproductive number will decrease. We ignore the 
discussion of mosquito’s susceptibility, thought it also increases with the basic reproductive 
number because it sound impracticable to think about vaccinating mosquitoes. However one 
can vaccinate mosquitoes and reduces the basic reproductive number. 

This study shows that reducing susceptibility of one population even up to zero does 


not guarantee the eradication of ZIKA disease. Figure 6(a) and 6(b) shows that the basic 





reproduction number is greater than one when either human beings susceptibility is zero or 
primate susceptibility is zero. This implies that when all human beings are either vaccinated 
or possess natural immunity against the disease, the disease will still circulating within the 
primates population, and having the possibility of reappear in human being population when 
significant number of not vaccinated individuals are born. 

Figure 7(a) and 7(b) display how the basic reproductive number varies with recovery rates. 
As the recovery rate of primate or human being increases, the basic reproductive number 
decreases. Since we made a mathematical assumption that recovery rate is a reciprocal of 
infectious period, we can reduce the basic reproductive number by reducing the number of 
infectious period, probably, through medicating the victims so as they can recover much 


earlier than before. 





Figure 7(c) shows the basic reproductive number decreases to zero asymptotically as the 
mosquitoes natural death rate increases. This suggests that we can also reduce the disease 
by killing large amount of Aedes mosquitoes responsible for ZIKA transmissions if ethics 
allow. Moreover, unlike vaccination and application of mosquito nets and condoms, the use 
of mosquito pesticides alone seems to be able to eradicate the disease if we can use them, 
heavily, to kill large amount of mosquitoes. By increasing mosquito’s birth rate from 1/14 
to 15/14, the basic reproduction number decreases from 4.1901 to 0.9571. This implies 


possibility of eradicating the disease by dealing with mosquitoes alone. 
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4. CONCLUSION AND RECOMMENDATIONS 


This work shows the role played by primates in dynamics of ZIKA disease. It reveals 
that it is not easy to reduce the basic reproductive number below one if we ignore the 
dynamics of ZIKA in primates population. It is observed that treating only human beings 
population cannot let the basic reproduction number become less than one. ‘Treating human 
beings population alone, may let the disease disappear in the population but still circulating 
in primates population and ultimately rebound back in the human beings population via 
mosquitoes. 

Our study generally, suggests that any control strategies - vaccination, the use of mosquito 
nets, pesticides, condoms, and abstinence from sexual activities - can be used to reduce ZIKA 
infection to some extent since they all have an impact on basic reproductive number. 

However, vaccinating human beings alone against ZIKA disease or the application of 
condoms, mosquito nets, and abstinence from sexual affairs alone can never eradicate the 
disease. Instead, the disease can be eradicated by application of the following possible 
(but not limited to ) combinations of different control strategies: vaccine both primates 
and human beings; heavily use of mosquito pesticides that can let mosquito’s death rate 
be at least 15times more than its natural death rate; a combination of the use of mosquito 
pesticides and the use of nets, condoms, and abstinence from sexual activities; a combination 
of the use of mosquito pesticide and vaccinating human beings/primates. 

This work does not suggest any best set of control strategies for fighting against ZIKA 


disease. It rather suggests some possible strategies that can be used to fight against the 








disease. We therefore, recommend an optimal control study, to be done, that will find out an 





optimal set of control strategies that can eradicate ZIKA disease effectively and efficiently. 
The study should specify exactly (but not limited to): a threshold proportion of human 
beings population or/and primate population to be vaccinated with/without combination of 
other strategies; a threshold mosquito birth rate to be attained with/without combinations 


of other strategies if viable. 
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